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Abstract 

We calculate the rates for the charged lepton flavour violating (LFV) decays 
£i —7- Ijj, T —7- ijT, T —7- irj^'^ iJL~ —7- e~e^e~, the six three-body leptonic de- 
cays —7- ij^^j'i'i^ and the rate for fi — e conversion in nuclei in the Standard 
Model (SM3) extended by a fourth generation of quarks and leptons (SM4), as- 
suming that neutrinos arc Dirac particles. We also calculate branching ratios for 
Kl — ^ fie, Kl — ^ TT^jie, B^^g l^e, B^^s and B^^g t/j,. We find that the 

pattern of the LFV branching ratios in the SM4 differs significantly from the one 
encountered in the MSSM, allowing to distinguish these two models with the help 
of LFV processes in a transparent manner. Also differences with respect to the Lit- 
tlest Higgs model with T-parity are found. Most importantly the branching ratios 
for Ci — >■ ijj, T — >■ in, T — >■ ir]^'\ ii~ — >■ e~e^e~ , t~ — )• e~e^e~ , t~ ii~ ji^ , 
T~ —7- e~ fL'^ and r~ — t- fi~e^e~ can all still be as large as the present experi- 
mental upper bounds but not necessarily simultaneously. Also the rate for ji — e 
conversion in nuclei can reach the corresponding upper bound. 



1 Introduction 



It is well known that in the Standard Model with three generations of leptons (SMS), 
the FCNC processes in the charged lepton sector, like ii — )■ £^7 and /x^ — )■ e~e+e~, are 
very strongly suppressed due to tiny neutrino masses. In particular, the branching ratio 
for /i — )■ 67 in the SMS amounts to at most 10~^^, to be compared with the present 
experimental upper bound 1.2 ■ IQ-^^ [l], and with the one that should be available in 
the coming years, ~ 10~^^ I^H^- Therefore any observation of lepton flavour violation 
(LFV) in the foreseeable future would be an unambiguous signal of New Physics (NP) 
beyond the SMS. 

Indeed, in extensions of the SMS like the MSSM, the Littlest Higgs model with T- 
parity (LHT) and Randall-Sundrum (RS) models, the presence of heavy leptons allows 
for much higher branching ratios for LFV processes, provided the mixing of these new 
leptons with the ordinary leptons is sufficiently large. On the other hand, none of these 
processes has been observed in experiment so far, and therefore we have to conclude that 
LFV rates are strongly suppressed in nature. As a consequence, these observations often 
put severe constraints on the parameter space of NP models in which heavy leptons are 
present. Moreover, these decays, being unaffected by hadronic uncertainties, allow for 
a clear distinction between different NP scenarios, in particular when several branching 
ratios are considered simultaneously and patterns of LFV in said extensions are compared 
to each other. For instance, a clear distinction between the MSSM and the LHT model [i] 
should be possible in this manner, once the relevant experimental data will be available. 

Models like the MSSM, LHT and RS scenarios contain many free parameters, and 
consequently do not always allow for clear-cut conclusions. On the other hand, one of the 
simplest extensions of the SMS, which contains heavy leptons and which has much less 
free parameters than the MSSM, LHT and RS scenarios, is the addition of a sequential 
fourth generation (4G) of quarks and leptons (SM4). 

In two recent papers [5|[6], we have analysed FCNC processes in the quark sector 
within the SM4, finding often sizeable departures from the SMS predictions for a multi- 
tude of rare decays of mesons and for — K^, B^^ — B^^^ — particle-antiparticle 
mixings, including in particular CP-violating observables. For recent work on SM4, see, 
for instance, [7 15 and references in (5|[6]. 

The goal of the present paper is a new analysis of the £j — )■ £^7, li — )■ Mk decays and 
of other LFV processes within the SM4, with the aim of finding the pattern of LFV in 
this model and constraining the masses of the new leptons and their mixing with the 
ordinary SMS leptons. There have been several analyses of LFV within the SM4 in the 
past 16 -18 , but to our knowledge, no detailed analysis of correlations between several 
LFV branching ratios has been presented to date. As we will demonstrate below, these 
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correlations could provide us with a very valuable tool for distinguishing the SM4 from 
other NP scenarios. 

Our paper is organised as follows. In Section [2] we provide the 4x4 leptonic mixing 
matrix Usma and discuss its various properties, taking into account the present experi- 
mental information about the 3x3 PMNS matrix. The subsequent section can be con- 
sidered as a compendium of formulae for the most interesting branching ratios for LFV 



processes within the SM4. In Section 3.1, formulae for the dipole transitions /i — )■ 67, 
r — 7- /i7 and r — )■ 67 are presented. The corresponding formulae for the three types of 
four-lepton transitions £~ — )■ ^J^'j^J, t ^7 ^7^7 r~ — )■ i'i^i' are presented in 



Section |3.2[ In Section |3.3| we turn our attention to the semi-leptonic r-decays, r — )■ fxP 
with P = TT, 77,77', and in Section 3.4 we calculate the fi — e conversion rate in nuclei. 



The Sections 3.5 and 3.6 deal with the decays K — )■ {TT^)fie and B^^s 2£, respec- 



tively, while in Section 3.7 the issue of (g — 2)^ is discussed. In Section [i] we present a 
detailed numerical analysis of all processes listed above, paying particular attention to 



several correlations. As demonstrated in Section 4.7, these correlations allow for a clear 



distinction of the SM4 from the MSSM and the LHT model. Finally, in Section |5] we 
conclude our paper with a list of messages from our analysis and with a brief outlook. 
Few technical details are relegated to the appendix. 



2 The 4x4 Leptonic Mixing Matrix t/sM4 

The most general leptonic mixing matrix in the 4G case contains 

6 mixing angles , 3 Dirac phases , 3 Majorana phases. 

A standard parametrisation is obtained by treating the mixing angles and Dirac phases 
in analogy to the quark sector (see e.g. [5, 19 ; for early work see 20-22]), with the 
Majorana phases contained in an additional diagonal matrix. 



UsMi = U34 14(524) f^24 h{-52i) hiSu) Uli /4(-5i4) U23 h{5i3) Ui3 hi-Su) Uu h 



Maj. 



(2.1) 

where Uij are rotations in the i-j plane, parameterised by corresponding mixing angles 
6ij, and 



[h{a)],k = e*°'- , /Maj. = diag[e*"S e^"^ e*"^ 1] 



(2.2) 



contain the Dirac and Majorana phases, respectively. The observables we are going to 
consider in the following are insensitive to the Majorana phases a,, which will therefore 
be dropped in the following. 
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Concerning the mixing angles and Dirac phases, it is well known that the SMS lepton 
sector behaves very differently as compared to the SMS quark sector. In particular, the 
PMNS matrix for SMS leptons is known to follow an approximate "tri-bi-maximal" 
mixing pattern [23], where 
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(2.3) 



which corresponds to the situation where the SG mixing angles satisfy S23 = l/"\/2, 

S12 = l/VS, Si3 = 0. 

A priori, it is not clear whether such a pattern could or should be extended to an 
SM4 lepton-mixing matrix, leading to potentially large mixing angles between the new 
lepton generation and the SMS ones. However, as we will see in more detail below, the 
current experimental situation already excludes large new mixing angles 6'j4 with the 4G 
leptons, and therefore we should rather consider 
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(2.4) 



as a starting point The deviations from this mixing pattern can then be conveniently 
described in terms of an almost diagonal mixing matrix. 



UsuA — 



residual 
SM4 



rrmax rMaj. 
'^SM4 ■ -' 



(2.5) 



where is parameterised in terms of small mixing angles Aij and S Dirac phases, 

and can be treated in an analogous way as the 4G quark mixing matrix. 

In particular, as in the quark sector, we may require the mixing angles Ajj to fulfill 
consistency relations [5 25 , 



AjfeAjfc < Aij 



(no summation over k) . 



(2.6) 



For k = 4, the product Aj4A 



i4 



lUi^U*^ on the left-hand side of this relation deter- 



mines, for instance, the size of radiative 



decays (see next section), which in turn 



^This ansatz reflects the special role of the fourth-generation neutrino, which requires some particular 
theoretical framework to be realized (see e.g. [24l). 
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set a lower bound (order-of-magnitude-wise) on the deviations of the PMNS matrix from 
tri-bi-maximal mixing, with 

While the deviations from tri-bi-maximal mixing Ajj - with the present experimental 
bounds - can still be of order 10 — 20%, the radiative LFV decays (see below) constrain 
the products 

A14A24 ^ 3.5 X 10-^ 
A14A34 < 1.8 X 10~\ 

A24A34 < 1.8 X 10-2. (2.8) 



3 Compendium for LFV in the SM4 



3.1 Dipole Transitions 

The diagrams for /i — )• 67 in the SMS and SM4 are shown in Fig. [T| and analogous 
diagrams exist for r — )■ /X7 and r — )■ 67. As demonstrated in |4], the relevant branching 
ratios can be found by inspecting the analogous calculation for B — )■ Xs'y. Neglecting 
tiny contributions from ordinary neutrinos and using formulae (3.8), (3.20) and (3.21) 
in we easily find 

2 
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where we have defined yj = ml./M^ and 
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The involved branching ratios of leptonic r decays are [26], 

Br(r" ^ UrC-Ue) = (17.84 ± 0.05)% , 
Br(r" ^ i^r/i"^^/.) = (17.36 ± 0.05)% , 

and the loop function 

H(x) = D'Jx) - Ie'Jx) , 



(3.9) 
(3.10) 
(3.11) 

(3.12) 

(3.13) 
(3.14) 

(3.15) 




is given in terms of functions D'q and E'q, known from the analysis of the B — )■ Xg'y decay, 
see Appendix [A} 

In writing (3.9)-(3.11 ), we have neglected electroweak (EW) corrections and non- 
unitarity corrections in the PMNS matrix that affect the leptonic decays and have been 
recently discussed by Lacker and Menzel [16]. Similarly, we have neglected corrections 
~ 0(m^/m^), ~ 0{m1/m^) and ~ 0{m?^/m'^^. These corrections amount to at most a 
few percent and would only be necessary in the presence of very accurate experimental 
branching ratios. 

The important virtue of formulae (3.9)-(3.11) is that, taken together, they allow for 
a direct determination of the ratios of the elements \Uei\, \U^i\ and |?7r4|, independently 
of the mass my^. In particular, we have 
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' Br(r~ yr^-'j^ji) 
Br( 



m4 



U, 



e4 



u, 



Br(r — > u^e z/g) . 



(3.16) 
(3.17) 
(3.18) 



For precision measurements, also the higher order corrections have to be included. 
Finally, let us note that in contrast to |16|, where the formulae of Altarelli et al. [27] 



were used, we are able to give the results (3.9)-(3.11) in an analytic form without any 
phase space integrals to be evaluated. 



3.2 Four-Lepton Transitions 

3.2.1 The Decays — >■ e^e+e", t~ — >■ fi~fj,'^fi~ and r~ — >• e~e'^e~ 

Next, we will consider decays with three leptons in the final state. Dipole operators, 
photon penguins, Z-penguins and box diagrams contribute here. In general, this will 
generate a rather non-trivial Dalitz distribution for the final state 28 which should 
be taken into account in the experimental analysis. In the SM4 analysis, we find that 
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the dipole operators are generally sub-leading, and therefore the Dalitz distributions are 
rather fiat functions of the invariant masses of the final-state lepton pairs. It is therefore 
sufficient to consider the formulae for the (partially) integrated branching ratios, which 
can be directly obtained from the corresponding expressions in [4] by appropriately 
replacing the loop functions by those of the SM4. We first find, using (5.8) of |4], 

r(/i^ —7- e^e+e^) 



Br(/i — )• e e~^e~ 



a 



(lie) 

XI 



3Z' + 3 ZH{y^) + H{y^y [log 
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+ 



sin^ 9i 



-ZY, 



2 sin 9w sin uw 
The one-loop functions entering this formula are easily found: 

K = Yo{y,)-\U,,\^So{y,), 



sin^ 9w 



-H{y,)Y^ 



(3.19) 

(3.20) 
(3.21) 

(3.22) 

The expressions for Yq, Sq, Cq, -Dq and Eq can be found in Appendix lAl and H{x) is 



Z = Co(i/4) + -^GQ{y,) 



where 



Gq{x) = Dq(x) 



'Eq{x) 



given in (3.15). Note that Dq(x) and Eq(x) diverge logarithmically for x — )■ 0, but Gq{x) 



vanishes for x — )■ 0. Since also Co(0) = 0, the contributions of light neutrinos can be 
neglected. 



For r — yU^/i , we make the following replacements in (3.19) 



Y.-^Y, 



(3.23) 



appropriately adjust the charged lepton masses, and multiply (3.19) by Br (r — > z/^/i 



The function Y^ can be obtained from Ye in (3.20) by replacing |t/e4p I^^mP- -'■^ 



case of r — > e e~^e , the corresponding replacements with respect to (3.19) are 



X4 



^Xa 



(3.24) 



and again the appropriate charged lepton masses have to be used. Now Br(r Ue) 
enters the formula for Br(r^ — )■ e~e^e~) as an overall factor. 



3.2.2 The Decays r — )■ e /i+e and r — )■ e+/^ 

Again following [4| and adjusting the formulae given there to the SM4, we find 



Br(r — )■ e /i^e ) 



1927r3 



47r2 



(3.25) 
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where summation goes over all four neutrinos. Using unitarity of the 4G leptonic mixing 

s to 

(3.26) 



matrix and the fact that only the 2,j = 4 contribution is relevant, (3.25) simplifies to 

2 



Br(r — 7- e ^ e 



1927r3 

where we have introduced the function 



(re) (lie) 



R2{x) = R{x, 0, 0, x) = 1 - 2P(x, 0) + P{x, x) 



(3.27) 



The branching ratio for r — )■ /i e'^jj, is obtained from (3.25 ) by interchanging n e 



and using x^j"^^ 



Xj 



3.2.3 The Decays r — >^ e+e and r e n'^fi 

Adapting the formulae of |4] to the SM4, we obtain 

with the differential decay rate K^^^s) 

„2 



(3.28) 
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The relevant Wilson coefficients read 



^r/i _ eA4 
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(3.30) 
(3.31) 



7W Siii L7pi/ 

where we have introduced the contribution due to additional box diagrams 



(3.32) 

In the above we neglected RGE running of a as well as operator mixing. The integral 



4 sin^Ow 



in (3.28) can then be performed analytically, and we arrive at 

„2 



Br(r — 7- fi e+e ) a 



-8\C^''mi-z) {8-z-z^) 
+2ARe{C^''C^''*) (1 - z)' - 48 IC^^'^p log (z) 



(3.33) 
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where we have introduced 



z = 



(3.34) 



Again, the formula for Br(r — t- e yU^/i ) is obtained from its analogous counterpart by 



interchanging ^ e in (3.33) 



3.3 Semi-Leptonic r-Decays 

The upper limits from Belle and Babar for the branching ratios of the decays r — )■ /ivr, 
r — 7- /ii] and r — )■ ^7]' are given in Table [2] Analytic expressions for these can be obtained 
directly from the corresponding formulae (4.12) and (4.17) in |4] by replacing the loop 
functions of the LHT model by the loop functions of the SM4. Neglecting suppressed 
pion and muon mass contributions of order (9(m^/m^) and (9(m^/m^), we find 



Br(r — )■ jin) 



X4 



(X + YY 



(3.35) 



1287r3 sin^ Ow 

with Tr and rrir being the lifetime and mass of the decaying r lepton. The branching 



ratio for the r — t- evr decay can be obtained very easily from (3.35) by simply replacing 
with f/e4- The generalisation of (3.35) to the decays r ^ fir] and r — )■ fir]' is slightly 
complicated by mixing in the r] — r]' mesonsystem and has been discussed in detail in jij. 
Proceeding as there, one obtains 



Br(r — )• nr]) 



1287r3 sin^ Ow 



X4 



2 / COS 6^8 Fg 



(X + Y) 



^sin^o^(X-2F) 



(3.36) 



Br(r —7- fir]') 



GWF^mlr^ 
1287r3 sin^ Ow 



X4 



2 / sin 6*8 Fs 



V3 F^ 



{X + Y) + 



^cos^o^(^ 



2Y) 



(3.37) 



where the mixing is described in terms of octet and singlet decay constants Fg, Fq and 
two mixing angles 9s, Oq 29 32 . Numerical input values are collected in Table [l] The 
functions X and Y are given by 



X = Xo(y4) 

Y = lo(l/4) 



+ 



i=d,s,b,b' 



(3.38) 
(3.39) 



i=u,c,t,t' 



Here Xi = rn^jM"^, rri^jM^ are functions of the quark masses, and Vij are the elements 
of the V4G mixing matrix in the quark sector. The functions F^P- and F"'^ are known 
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from our analysis of rare K and B decays in |5]. We recall their explicit expressions in 
Appendix lAl 



3.4 ^ — e Conversion in Nuclei 



Here, as in ^ , we use the general formula (58) of 33 to find the approximate conversion 
rate 



r(/iX -> eX) 



(lie) 



(3.40) 



X 



Y 



Z{AZ + H{y,))-{2Z + N)^^ + {Z + 2N) . 

sm U\Y sm Uy/ 



where X and Y are given in (3.38) and (3.39), and H[x) and Z are given in (3.15) 



and (3.21). Z and N denote the proton and neutron number of the nucleus. has 
been determined in [34|]37| and F{q'^) is the nucleon form factor. For X = 22Ti, one has 
Zeff = 17.6 and F{q'^ 



0.54 



38 



The /i — e conversion rate R(/iX — )■ eX) is then given by 

r(/iX ^ eX) 



R(/xX ^ eX) 



capture 



(3.41) 



with r^p^y^g being the /x capture rate of the element X. For titanium the experimental 



value is given by 39 



pTi 

capture 



(2.590 ± 0.012) ■ lOS"^ (3.42) 

In our numerical analysis of Section |4] we will restrict ourselves to /i — e conversion 
in 22^1, for which the most stringent experimental upper bound exists and where the 
approximations entering ( |3.40[ ) work very well. For details, we refer the reader to the 
discussion presented in [33|[38|[40] . 



3.5 Kl^s and Kl^s n^fie 

In the SM3 the decay Kl — ?■ yue can proceed through box diagrams in the case of non- 
degenerate neutrino masses, but similarly to /i — )■ 67 its rate is too small to be measured. 

Also in the SM4, Kl — )■ fie proceeds through box diagrams as shown in Fig. [2| but 
due to the large mass of the 4G neutrino this contribution now becomes relevant. The 
effective Hamiltonian corresponding to these diagrams is given by 

Hes = VAf )xr^i'(^-^^)(^^)v-A(eyu)y^A, (3.43) 
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Figure 2: Diagrams contributing to the Kl — )■ /xe decay in the SM4. 



where 



(3.44) 

and the function P{xi, yj) can be found in Appendix [A} 

Using the unitarity of the 4G leptonic mixing matrix and the fact that i/j ~ for 
j 7^ 4, this expression simphfies to 

Heff = ^MlxT^ Y. ^S""^ [^(^- - P{xi, 0)] {sd)v-Amv-A ■ (3.45) 



The evaluation of Bt{Kl — )■ yue) by means of Hes in (3.45) proceeds along the deriva- 
tion presented in Section 7 of 13. We accordingly find 



(lie) 



X 



Re{Xf^)R{x„0,0,y,) 



\i=c,t,t' 

where the function R is defined as follows 



(3.46) 



R{xi, Xj, yk, yi) = P{xu yi) + P{xj, yt) - P{xj, yi) - P{xu yt) 



and 26,41 42 



Br(fs:+ ^ /i+z/) = (63.44 ± 0.14)% , 



(3.47) 



= 4.117 ±0.019, iKsl = 0.225 ±0.001. 

(3.48) 

Note that in general \Vus\ from a SM3 fit of semileptonic K decays is not longer valid 
in the SM4. However since a reanalysis of this fit in the context of the SM4 is clearly 
beyond the scope of the work, we use the above value for simplicity. Experimentally we 
have |43l|44l 



Bt{Kl fie) = Bt{Kl /i+e") ± Bt{Kl /i"e+) < 4.7 ■ 10 



-12 



(3.49) 
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Similarly, following 14] we find 



Bt{Kl 71^ fie) = Bt{Kl vrV+e-) + Bt{Kl TrV^e^ 



r{K+) IK, 



(ue) 

xT 



(3.50) 



\i=c,t,t' 



where |2& 



Br(K+ ^ ny+iy) = (3.32 ± 0.06)% . 



(3.51) 



3.6 Lepton-Flavour Violating B Decays 



A detailed study of the decays B^^s /^e, B^^s — )■ re and B^^s — )■ r/i in the LHT model 
has been presented in Section 8 of [4] . The formulae given therein can easily be adapted 
to the SM4 case and are summarised by the following two equations: 



Br(5, ^ 



GlM^ t{B,) 



Br(5+^£+z/,,; 



Br(5, ^ hh) 



167r^\Vub\^ t{B 

J2 Af)i?(x„0,0,i/4) 

■.=c,t,t' 



A4 



(3.52) 



167r4|KbPr(i?+)MB,F. 



2 Br(S+^£+i/,,; 



X4 



(^1^2) 



5^ X'f^R{xi,0,0,y,) 



i=c,t,t' 



(3.53) 



where £i and £2 denote the two leptons in the final state with m^^ > m^j. For our 
numerical analysis we used the SM predictions for Br(i?+ — fi^iy^) and Br(i?+ — r^Ur) 



45,46 



Br(fi+ ^ /i+z/^) = (3.8 ± 1.1) ■ 10' 



(3.54) 
(3.55) 



Br(5+ T+Ur) = (0.8 ± 0.12) ■ 10"^ . 
This is necessary, because currently only Br(i?+ — )■ r^Ur) = (1.67 ±0.39) • 10^"^ is 



measured 47-49 . The SM3 value is on the lower side of this measurement but still 



consistent within errors. 
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3.7 Anomalous Magnetic Moment of the Muon 

The one loop contribution to = {g — 2)^/2 = (0^)3^43 + {ci^l)sMA SM4 can be 

obtained in analogy to our derivation in . Only the triangle diagram with a W boson 
and two heavy neutrinos running in the loop is relevant here. Adapting the corresponding 
expression (11.11) in fi] to the SM4 by dividing by a factor f^/(4/^), removing the sum 
over different flavours and adjusting the matrix elements jV^'^p — )■ |f/;i4p we obtain 

MsM4 = - ^^2^ ^l I^^mI^ (^2(1/4) - L2{0)) , (3.56) 

where the function L2{x) is given in Appendix |Aj Two comments are in order at this 
point: 

1. Since -^2(2/4) — i^2(0) > 0, the SM4 contribution tends to decrease and thus 
pushes it even further away from the experimental value. 

2. It turns out that after imposing the constraints from lepton universality and ra- 
diative decays, the SM4 contribution to becomes negligible compared to the 
theoretical uncertainties. 



4 Numerical Analysis 



4.1 Preliminaries 



rrip 



0.5110 MeV 



= 105.66 MeV 
= 1.77684(17) GeV 
= 290.6(1.0) ■ lO-^'ps 

Mw = 80.425(38) GeV 

a = 1/137 

Gf = 1.16637(1) • 10-5 GeV 
sin^^VF = 0.23122(15) (26 



-2 



192.8(9.9) MeV 
238.8(9.5) MeV 



50 



r(5,)/r(fi+) = 0.934(7) 
t{Bs) = 1.425(41) ps 
t{B+) = 1.638(11) ps 
Mb^ = 5.2794(5) GeV 
Mb, = 5.3675(18) GeV 

\Vub\ = 3.68(14) ■ 10-3 (iel 



Fs/F^ = 1.28 (ChPT) 
Fo/F^ = 1.18(4) 
^8 = -22.2(1.8)° 
60 = -8.7(2.1)° 
F^ = 130 ± 5 MeV 



51 



26 



Table 1: Values of the experimental and theoretical quantities used as input parameters. 



The great simplicity of the analysis of LFV within the SM4 when compared to NP 
scenarios such as the general MSSM, the LHT and RS models is the paucity of free 
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parameters. The analysis also simplifies compared to the quark sector since the contri- 
butions of SMS leptons in loops can be neglected, except when they are relevant in the 
context of the GIM mechanism. 



We note that 
quantities 



(lie) 



j7 and the decays to three leptons are fully governed by the 



(re) 

x\ 



Xl 



e4| ) 



If/, 



(4.57) 



ments 26 



> 90.3 GeV 



and calculable functions of the neutrino mass m^,^ which is bounded by direct measure- 

(95% C.L.) . (4.58) 

Therefore strong correlations between the ti — )■ and — )■ Mk decays are to be ex- 
pected. While this expectation will be confirmed in the course of our numerical analysis, 
we will see that the possible ranges for various observables entering these correlations 
will still be rather large. 

Semileptonic decays and fi — e conversion in nuclei involve also parameters in the 
quark sector that enter through box diagram contributions to the functions X and Y 
in (3.39) and (3.38). These contributions are however constrained through our analysis 
of the quark sector in the SM4. In our analysis we also take into account constraints 



present outside the LFV sector, in particular those from 16 



For our numerical analysis of processes involving quarks we used the points of our pre- 
vious analysis Our parameter points were generated using uniform random numbers, 
and we explicitly do not assign any statistical meaning to the point densities. We in- 
cluded the effect of a modified Fermi constant Gp due to the breaking of three-generation 
lepton-universality 16 and included the decays r — t- /xz/^z/t- and r — )■ ez/gZ/T- to constrain 
the parameters. Contrary to 16 we do not find a significant effect of the K31 decays, 
but this is due to our much more conservative error treatment. On this note we want to 
reemphasise the need for a consistent fit of the EWP data, CKM matrix elements from 



semileptonic decays, Gp and similar well known inputs 15,16,19,52 54 , in the context 
of the SM4. 



4.2 fi — )- e 7, /X —^e e+e and — e Conversion 

In Fig. [3] we show the correlation between yU — )■ 67 and yU~ — )■ e~e^e~ together with the 
experimental bounds on these decays. We observe: 

• Both branching ratios can easily reach the present experimental bounds in a cor- 
related manner. 

• However, for a fixed value of either branching ratio, the other one can still vary 
over one order of magnitude. 
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Br(;j ^ e y) 



Figure 3: Correlation between Br(/i — )■ 67) and Br(/i — )■ e e+e ). Points that agree 
with the currently measured fi — e conversion rate in 22Ti (3.42) are shown in blue, while 



gray points violate this bound. The shaded area indicates the projected experimental 
bound on Br(/i — )■ 67) from the MEG experiment at PSI. 



Next in Fig. |4] we show the correlation between the — e conversion rate in 22Ti and 
Br(/i — )■ 67), after imposing the existing constraints on Br(yU — )■ 67) and Br(/i~ — )■ e~e"*'e~). 
We observe that this correlation is weaker than the one in Fig. [3] as now also quark pa- 
rameters enter the game. Still, for a given Br(/i — )■ 67) a sharp upper bound on the 
/i — > e conversion rate is identified. Furthermore, we find that the /i — e conversion rate 
in titanium is generally larger than the current experimental bound, but the bounds 
on both branching ratios can be simultaneously satisfied. Yet it is evident from this 
plot that lowering the upper bounds on both observables in the future will significantly 
reduce the allowed regions of the leptonic parameter space in the SM4. 



As pointed out by 16 , the combination of results from leptonic r decays and radiative 
H decays efficiently constrains the involved PMNS parameters \Ue4\ and \Ufj_4\. We show 
these bounds in Fig. [5| adding the constraint from n — e conversion which turns out to 
be the most stringent one. 

4.3 The Decays r — >^ and r — )■ 67 



In Fig. |6]we show the correlation between Br(r — fi'y) and Br(r — 67), imposing the 
experimental bounds on /i — )■ 67 and /i~ — > e^e+e". We observe that they both can be 
individually as high as few times 10~^ and thus in the ball park of present experimental 
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Figure 4: Correlation between Br(/i — )■ 67) and R(/iTi — )■ eTi). The shaded areas 
indicate the expected future experimental bounds on both observables. 




Figure 5: Constraints on the allowed range of |f/e4| and |t//i4| resulting from lepton 
universality {la/ 2cr/3cr: dark gray/gray/light gray area, respectively) and the current 
experimental bounds on /i — )■ eee, /i — )■ 67, and /x — e conversion (thick black lines). 
The contour lines indicate the ratio Gp^^'^/Gp'^^, where Gp'^'^ is the value of the Fermi 
constant extracted from muon lifetime measurement assuming 4 generations, and Gp^^ 
is the usual SMS Fermi constant (s. [Ig]). 
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Br(T ^ e y) 



Figure 6: Correlation between Br(r — )■ fi'j) and Br(r — 67). 



upper bounds. The maximal values however cannot be reached simultaneously, due to 
the /i — 67 constraint. Thus finding both branching ratios at the 10~^ — 10~^ level will 
basically eliminate the SM4 scenario. 



4.4 The Decays r fiTZ, firj, firj' and r — >^ fi'j 

In Fig. [7] we show Br(r — t- fin) a function of Br(r — )■ /i7), imposing the constraints from 
fi ^ e'j and fi^ e^e^e^. We find that Br(r — )■ /ivr) can reach values as high as the 
present experimental bounds from Belle and BaBar, which is in the ball park of 10^^. It 
is evident from (3.35)-(3.37) that Br(r — )■ firj') and Br(r — )■ fir]) are strongly correlated 
with Br(r — ?■ yuvr), so we choose not to show the respective plots for these processes. 

Completely analogous correlations can be found also for the corresponding decays 
r — )■ evr, er], erj and r 67. Indeed, this symmetry between the t ^ fi and r — i- e 
systems turns out to be a general feature of the SM4, that can be found in all decays 



considered in the present paper. We will return to this issue in Section 4.7 



An immediate consequence of these correlations is that the observation of a large 
r — 7- /i7 rate will immediately imply a large r — )■ /ivr rate and vice versa. Still, for a fixed 
value of either branching ratio the second one can vary by almost an order of magnitude. 
Analogous statements apply to r — )■ fi{e)ri and r — t- /i(e)?7'. 
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Br(r y) 



Figure 7: Br(r — )■ fin) as a function of Br(r — )■ fi'-f) 
4.5 Xx, — )- //e and Kl — >- 7r%te 

In Figs, [s] and [9] we show the results for Bt{Kl — )■ /le) and Bt{Kl — t- Ti^jie) as functions 
of Br(/i — )■ 67). Again strong correlations between these branching ratios are observed 
but the maximal values for Bt{Kl fJ'^) and Bt{Kl — )■ n^fie) are by several orders of 
magnitude below the present experimental bounds. 




Br(p ^ e y) 



Figure 8: Bt{Kl — )■ /le) as a function of Br(/i — 67). Colour coding defined in Fig. [3| 
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10-14 
Br(/j 



10-13 



10-12 



10 



-11 



Figure 9: Bi{Kl — t- 7r°/ie) as a function of Br(yU — )■ 67). Colour coding defined in Fig. ^ 



4.6 Upper Bounds 

In Table [2] we show the maximal values obtainable in the SM4 for all branching ratios 
considered in the present paper, together with the corresponding experimental bounds. 
We observe: 

• The branching ratios for eleven of decays in this table can still come close to the 
respective experimental bounds and as we have seen in the previous plots they are 
correlated with each other. 

• The remaining branching ratios are by several orders of magnitude below the 
present experimental bounds and if the SM4 is the whole story, these decays will 
not be seen in the foreseeable future. 

• Comparing to the results obtained in the LHT model for a NP scale / = 1 TeV 
[4,55, the SM4 allows for much larger branching ratios but the difference is much 
smaller for / = 500 GeV. 

We have also investigated the effect of additionally imposing R(/iTi — )■ eTi) < 5 ■ lO-^^ 
as a constraint, which we have chosen slightly above the experimental value 4.3 ■ lO-^^ 
in order to account for the involved theoretical uncertainties. We find that all maximal 
values collected in Table |2] depend only weakly on that constraint. 
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maximal value 


exp. upper bound 


/i — )■ 67 


1.2 ■ 


10"" (6.8 


■ 10^ 


-12) 


1.2 


■ 10 


1 1 
— 11 


H 


/X — )■ e e^e 


1.0 


1 r\ 19 / -1 

■ 10 (1 ■ 


10- 


12) 


1.0 ■ 


10- 


12 


|56] 


K(/uii — > eii) 


6.6 


■ 10-" (5- 
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4.3 ■ 


10" 


IZ 


1571 


r — 7- 67 
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■ 10-8 (3.9 


• 10- 


Q \ 
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10" 


Q 

-o 
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Q 
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10" 


Q 
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Table 2: Maximal values for LFV decay branching ratios in the SM4, after 
imposing the constraints on Br(/i — t- e7) and Br(yU- — t- e"e"'"e"). The num- 
bers given in brackets are obtained after imposing the additional constraint 
R(/iTi — )■ eTi) < 5 ■ 10"i^. The current experimental upper bounds are also 
given. 



This finding justifies that we did not take into account this bound in our numerical 
analysis so far, as it has only a minor impact on the discussed observables. We would like 
to stress that the maximal values in Table [2] should only be considered as rough upper 
bounds. They have been obtained from scattering over the allowed parameter space of 
the model. In particular, no confidence level can be assigned to them. The same applies 
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to the ranges given in Table [3] for the SM4 and the LHT model. 



4.7 Patterns of Correlations and Comparison with the MSSM 
and the LHT 



In [4 ,55 a number of correlations have been identified that allow to distinguish the LHT 
model from the MSSM. These results are recalled in Table [H In the last column of this 
table we also show the results obtained in the SM4. We observe: 

• For most of the ratios considered here the values found in the SM4 are significantly 
larger than in the LHT and by one to two orders of magnitude larger than in the 
MSSM. 

• In the case of — )■ e conversion the predictions of the SM4 and the LHT model 
are very uncertain but finding said ratio to be of order one would favour the SM4 
and the LHT model over the MSSM. 

• Similarly, in the case of several ratios considered in this table, finding them to be 
of order one will choose the SM4 as a clear winner in this competition. 



ratio 



LHT 



MSSM (dipole) 



MSSM (Higgs) 



SM4 



Br(/^ — >e e+e ) 
Br(/j— ^67) 

Br(r~— ^-e^e+e") 
Br(r-s>e7) 

Br(T— >/i7) 



0.02. ..1 
0.04. . .0.4 
0.04. . .0.4 



6 • 10-3 
1 ■ 10-2 
2 ■ 10-3 



~ 6 ■ 10-3 
~ 1 ■ 10-2 
0.06... 0.1 



0.06. 
0.07. 
0.06. 



2.2 
2.2 
2.2 



Br(T — >e IJ,'^ fJ. ) 



Br( 



>e7) 



Br(r e+e ) 

Br(r^— >e^e+e^) 
Br(T~— >e~/j+/^~) 

Br(r~— ) 
Br(-; — — >/j~e+e^) 



0.04. ..0.3 
0.04. ..0.3 
0.8. ..2 
0.7. ..1.6 



2 • 10-3 
1 • 10-2 
~ 5 
~ 0.2 



0.02... 0.04 
~ 1 ■ 10-2 
0.3... 0.5 
5... 10 



0.03. 
0.04. 



1.3 
1.4 



1.5... 2.3 
1.4... 1.7 



R(AtTi-5-eTi) 
Br(/i— >e7) 



10-3... 102 



5 • 10- 



0.08... 0.15 



10-12 ___26 



Table 3: Comparison of various ratios of branching ratios in the LHT model 55 , the 



MSSM without 63,64 and with significant Higgs contributions 65 ,66 and the SM4 
calculated here. 
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5 Conclusions 



In the present paper we have calculated branching ratios for a large number of charged 
lepton flavour violating decays in the Standard Model extended by a fourth genera- 
tion of quarks and leptons, assuming that neutrinos are Dirac particles and taking all 
presently available constraints into account. Our main messages from this analysis are 
the following: 

• The branching ratios for £i — ij'j, r — )■ in, r — t- £r]^'\ fi^ e^e+e^, — i- e^e+e^, 

— i- T~ — )■ e~^^^^ and r~ — > fi^e^e^ can all still be as large as the 

present experimental upper bounds but not necessarily simultaneously. 

• The correlations between various branching ratios should allow to test this model. 
This should be contrasted with the SMS where all these branching ratios are un- 
measurable. 

• The rate for n — e conversion in nuclei can also reach the corresponding upper 
bound. 

• The pattern of the LFV branching ratios in the SM4 differs significantly from the 
one encountered in the MSSM, allowing to distinguish these two models with the 
help of LFV processes in a transparent manner. The same statement applies to 
the LHT, as can be clearly seen from Table |3} 

• The branching ratios for Kl — )■ /ie, Kl — n^fie, Bd^s — )■ yue, Bd,s — ^ tc and 
Bd,s TyU turn out to be by several orders of magnitude smaller than the present 
experimental bounds. 
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A Relevant Functions 



In this appendix we collect the various functions entering the theoretical formulas for 
the LFV decays discussed in Sec. [3j 



Co{xi) 
Eo{xi) 



Xi Xi 6 3xi ~h 2 
Xi - 1 [x^ - ly 

4 , -19x3 + 25x2 x2(5x2 - 2xi - 6) , 

— logXj H -r^ + — — — — - logXi , 

9 ^ 36(xi-l)3 18(xi-l)4 ^ ' 

2, x,2(15- 16xi + 4x2) , Xi(18-llXi-x2) 

- 3 + 6(^7^1? '"^ + 12(l-x,.)3 



(3xf - 2xf ) 



2{Xi - 1)4 
3x? 



log Xi + 



ixf + 5xf — 7xj 
12(x, - 1)3 



2(Xi - 1)4 

X. 



(xf - 5x^ - 2xj 
logXi+ ^ ' 



4(x, - 1)= 



Yo{xi) 



Xi 



Zo{xi) 



Xi + 2 3xi - 6 

Xi - 1 (Xi-1)^ 

3xi 



Xi 4 



+ 



log Xi 



) [xi - 1 (Xi - 1)2 

1 , 18x4 - 163x3 + 259x2 - 108xi 
■9 + 144(x. - 1)3 

32xf - 38x3 _ 153,2 ^ ig^, 

log Xi . 



72(x, - 1)4 

For arbitrary arguments Xi,Xj, the function 5'o(xi,Xj) is given by 67 



^o('^25^j) XiXj 



(4 — 8xj + x2) log Xi 



+ j) - 



(A.l) 
(A.2) 
(A.3) 
(A.4) 
(A.5) 
(A.6) 
(A.7) 



(A.8) 



^ * 4(Xj - 1)2 (Xj - Xi) ' ^" 4(xi - l){Xj - 1) 

In the limit of e — )• in 5'o(xi + e, Xi — e) one recovers the SM3 version of 5*0 (xj), 



(A.9) 



Xi -4 + 15xi - (12 - 61ogXi)x2 + x3 

^0(X.) = - ^—y, . (A.10) 



The functions entering the Kl,s decays discussed in Sec. 3.5 read 



F^^(x„ 2/4) = B^'^ixi, 0) + 5^^(0, 2/4) - B^'^iO, 0) - E'^^lxi, 2/4) 

= -5'o(xi,?/4) , (A.ll) 
F-'^ix,, z,) = B'^ix,, 0) + B-^'iO, z,) - B-^'iO, 0) - B'^'ix,, z,) , (A.12) 
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with 



and 



U{Xx,X2) 
U{xi,X2] 



x\ log Xi 



+ 



^2 log X2 



+ 



1 



(xi - X2)(l - a;i)2 (x2 - - ^2)^ (1 - a;i)(l - X2) 
Xiloga;i X2logX2 1 



(xi - X2)(l - (X2 - Xi)(l - ^2)^ (1-Xi)(l-X2) 

We also encounter the function 

1 / . 7 \ xf log Xi 



(A.13) 
(A. 14) 

(A.15) 
(A. 16) 



(l-x,)(l-%) 
{xi - yj){l - yjY 



1 - 2xi + 



{xi - ?/j)(l - a;*)^ 



(A.17) 



Finally, the result for the anomalous magnetic moment of the muon is expressed in terms 
of 



Lo(x) 



6(1 -x)^ 



(-10 + 43x - 78x2 ^ 4g^3 _ 4^4 _ ^g^3 jQg^j 



(A.18) 
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